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1. INTRODUCTION 


Thermally sprayed coatings have been extensively used to enhance materials 
properties and provide surface protection against their working environments in a number 
of industrial applications. Thermal barrier coatings (TBC) are used to reduce the thermal 
conductivity of aerospace turbine blades and improve the turbine overall thermal efficiency. 
TBC allows higher gas operating temperatures and lower blade material temperatures due 
to the thermal insulation provided by these ceramic coatings. In the automotive industry, 
coatings are currently applied to a number of moving parts that are subjected to friction and 
wear inside the engine such as pistons, cylinder liners, valves and crankshafts to enhance 
their wear resistance and prolong their useful operation and lifetime. 

Recently, hexavalent chromium associated with hard chromium plating, was 
classified as a human carcinogen and environmental pollutant. Aerospace industry has 
traditionally used hard chrome plating as a corrosion protection material. Accordingly, 
NASA is developing thermal spray coatings to replace the hard chrome plating that is 
currently being utilized on the space shuttle main engine. These coatings are made of 
tungsten carbide and cobalt, chromium oxide, and FerroTic (iron based material with 
titanium carbide) to provide protection against wear, corrosion and hydrogen embrittlement 
of the low pressure liquid hydrogen carrying ducts on the shuttle main engine. One of these 
sections is the Low Pressure Fuel Turbo Pump (LPFTP) discharge duct used on the shuttle 
P^ain Thp^urt ramps lifljiujjwrirogenj^ of -253 0 C i-42 3 °F ) to the 


velocity of molten particles [2]. Critical review of the literature on thermal spray coatings, 
indicated that high velocity thermal spray processing conditions offer the best potential to 
minimize the occurrence of manufacturing related flaws and provide high quality coatings 


In recent years the High Velocity Oxygen Fuel (HVOF) system has been considered an 
asset to the family of thermal spraying processes Fig. (1). Especially, for spray materials with 
melting points below 3000 °K it has proven successful, since it shows economic advantages 
when compared to other coating processes that produce similar quality coatings [4], In such 
systems that produce high velocity particles, the oxide content of the coatings did not 
correlate to the particle temperature or the excess oxygen in the lean fuel conditions, but it 
is related to high substrate temperatures [5]. It is believed that the primary mechanism for 
the formation oxide inclusions occurs after a particle splat, when the hot coating is exposed 
to the oxygen contained in the relatively low velocity boundary layer. The exposure times 
of a given splat to the boundary layer are on the order of seconds, a factor of 10 3 longer than 
during particle flight [6]. The usual evaluation of optimized spraying parameters, which 
include stand off distance, fuel and oxygen flow rates, powder size, and barrel length, are 
based on numerous, extensive, and expensive experiments laid out by trial-and-error or 
statistical design of experiments and Taguchi methods [4], 
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Figure 1. Thermal Spray processes 


The application of simulation techniques to thermal spray processes has grown 
steadily over the past years because of the relatively inexpensive parametric analysis and 
operations. Once a model of the process is established and validated, system parameters 
become evident. The obtained fundamental understanding of the optimization process can 
be accomplished with a numerical computer model [4]. 

The objective of this investigation was to develop a computer model that simulates 
the thermal and gas dynamics of both gases and particles associated with the HVOF process 
to provide predictions of the particles velocities and temperatures. In the HVOF process 
oxygen and atomized kerosene are injected coaxially into the combustion chamber. Fig. (2), 
where they are mixed and ignited by a spark plug. The hot combustion gases are accelerated 
to supersonic conditions in a conversion diversion nozzle. At the exit of the nozzle the 
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powder is transversely injected with an argon or nitrogen earner gas. The powder is heated 
and accelerated in the gun barrel and jet region by the hot gas and eventually is deposited on 
the substrate. The model was partially validated against the very limited experimental data 
collected from Hobart Tafa using the JP-5000 gun and shown in Figure (2). 
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Fig. (2) Hobart Tafa HVOF/JP-5000 - with showen Numerical Values for Gas and Particle 
Velocity and Temperature 

2. MATHEMATICAL MODEL 

The current model development was accomplished in four main stages as described in the 
following : 

First Stage: Modelling of the gases flow inside the HVQF system without initial 

consideration to th e powder flow 171: 

The generalized flow equations with the influence parameters were used to model the gas 
flow in the HVOF system shown in Figure (2). The temperature and Mach Number equations 
(1). (2) were numerically integrated for a single phase non-adiabatic friction flow with 
variable specific heats [7], 
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thermal spray system for the inlet conditions given in Fig. (1). Due to the appearance of the 
term (1-M 2 ) in the denominator of equations (1) and (2), the numerical model experienced 
some degree of singularity in the vicinity of the throat section. The model quickly recovered 
as the solution propagated downstream of the throat area. However, this singularity problem 
could be completely rectified by the application of a Taylor series expansion on these terms 
in the throat proximity zone. 

Second Stage: Empirical Correlation for the mean axial velocity and temperature decay of 
a free supersonic , iet plume : 

Visual studies such as Schlieren flow visualization showed a "potential core" with 
an embedded shock-diamond structure that is formed in the supersonic zone due to the 
underexpanded nature of the jet, Fig (3), [6], Further down stream large turbulent eddies 
were generated by the large velocity and temperature gradients at the boundary layer between 
the jet and the ambient air. This mixing zone is further subjected to the effect of a large 
density difference between the hot jet core and the comparatively cold and slow ambient 
atmosphere. Small particles 5 pm or smaller will fully track the turbulent motion of the 
fluid, however, much larger particles, such as typical HVOF metal spray powders, are 
generally unaffected by the eddies and remain in the relatively high temperature, low density, 
and least motion resistance zone near the jet centerline. 
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Fig. (3) Supersonic Jet Structure 

In an attempt to better understand the noise generation mechanism in supersonic free 
jets, considerable investigation efforts have been devoted to the measurements of various 
flow parameters and the study of how these quantities vary with jet flow conditions [8]. 
Earlier mean velocity measurements [8] using laser velocimeter and hot-wire anemometer 
resulted in an empirical formula which gives the variation of the potential core length. Figure 
(3) with the Mach number for both heated and cooled jets, it read as follows: 
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Another general empirical correlation for the supersonic axial mean velocity decay with 
range of validity between Mach Number 0.3 - 1 .4., this correlation read as follows [8]: 
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In the HVOF applications the jet mean axial velocity and temperature ranges are between 
(1.5 - 2.5 M) and (1500 - 3000° K) respectively. The recent study of the heated jets indicated 
that the potential core length decreased as the jet temperature increased [9]. Thus, 
correlations (3), and (4) were modified to expand their range of validity by correlating recent 
measurements of the mean axial decay velocity and temperature obtained from NASA/ 
Langley [9] and University of Toronto [10]. Figures 4 and 4A show the mean velocity and 
temperatures correlations in relation to the measured values respectively. The comparison 
of these Figures 4, and 4A for free axisymmetrical jets indicated that the velocity decays 
faster that the temperature. 

Third Stage: Momentum Transfer Mechanism Between Gas and Particle 

The momentum equation, for either solid or liquid particles, are solved in a 
Lagrangian frame of reference moving with the particles. The equation of motion for the 
particle is written as [1 1]: 


"V TT = 7 p . A » c ° <V v p>I W - u p v p < 5 > 

at c. 


Where m p is the mass of the particle and V p is the velocity vector of the particle, C D is the 
drag coefficient, and the particle motion in two-phase flows depends upon the gas properties, 
particle properties, and p g , V g , P are the density, velocity and pressure of the gas, 
respectively. A p is the particle surface area and V p is the particle volume. All particles are 
assumed to be spherical. This equation of motion for a particle accounts for the 
acceleration/deceleration of the droplet, due to the combined effects of drag from the gas 
flow, and local pressure gradients in the gas. Because the gas flow pressure change is small 
in the friction flow inside the gun barrel and the free jet plume, the gas pressure gradient 
effect on the particle motion is neglected in comparison to the drag force. A literature survey 
in the area of two-phase flows reveals that a number of drag coefficient equations have been 
used to calculate particle motion in supersonic flows [12]. Because the particles are initially 
injected into a supersonic flow, each particle will have a shock wave on the upstream side. 
The drag induced on these particles in this supersonic flow is calculated with the following 
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Figure (4> Modified Axial Decay Velocity in Supersonic Jet Plumes 
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Figure ( 4 ^ Modified Axial Decay Temperature For Supersonic Jet Plumes 
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empirical correlation: 
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The drag coefficient for the particle is based on the local Reynolds number of the particle and 
is evaluated as: 
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here p is the molecular viscosity of the gas. The variation of the gas viscosity was evaluated 
from the graph and empirical formula shown in Figure (5). The following correlations have 
been found to be valid for a wide range of Reynolds number [11]; it reads as follows: 
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C D = 0.44 for R e > 10 3 


Fourth Stage: Heat Transfer Mechanism Betwee n The Gas and The Powder; 

The coefficient of heat transfer (h) between the particle and the gas can be determined 
from the following Ranz-Marshall semiemperical equation: 
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Where c p is the specific heat of combustion gases given by correlations shown in table (1), 
dp is particle diameter, and k g is the coefficient of gases thermal conductivity - it was 
evaluated from a correlation as a function of the temperature as shown in Figure (6). 
Considering that the particle maintained its spherical configuration all through the process, 
the temperature of the particle was calculated from the following expression: 
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3. RESULTS AND ANALYSIS 


The model predictions showed an expected sharp increase in both of the particle 
velocity and temperature when the powder was initially injected in the barrel, Figures (7A) 
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Fig.(7A) Predicted Gas and Particles Velocities Inside The JP-5000 Barrel 



Fig. (*}B) Predicted Gas and Particles Velocities In The Free Jet Plume of a HVOF System 
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Fig.(8A) Predicted Gas and Powder Temperatures inside The JP-5000 Gun Barrel 



Fig. (8B) Predicted Gas and Powder Temperatures in The Jet Plume 
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6. RECOMMENDATIONS 


(1) Measurements of velocity and temperature of gas and particle are necessary for the 
model validation 

(2) The model should be used to perform parametric study and system optimization 

(3) Rectify the singularity problem due to the appearance of the term (1-M 2 ) in the 
denominator - Use Taylor series expansion 
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